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- ABSTRACT 


Experimental studies were conducted on unfilled visco- 
mits re Materials to evaluate the fracture postulate which 
merece that crack vellocity is a unique function of stress 
i@reensity factor. HOnComtEnis, fracture characterization 
meses were conducted for constant load, constant displace- 
Tome, eanad COnstant displacement rate histories. A cyclic 
Stuay, a more sophisticated application of the constant dis- 
meteement rate Case, was also conducted as a further evalua- 
mom OL the £L£racture postulate. 

tiemaidinatke al theory used in obtaining viscoelastic 
stress intensity factors was based upon linear theory and 
sharp crack geometries. Moderate correlation was obtained 
gmonuge all test cases, notwithstanding the presence of 


Begeiricant Specimen strains and finite crack tip radii. 
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’ IB INTRODUCTION 


Over the years, fracture studies have led to the devel- 
@oment Of relationships among geometry, loading and crack 
Jength that prevail at the incidence of rupture and the 
Semsea Or crack propagation. It has been found useful to 


memmme a Stress intensity factor, K which increases as 


1! 
mmemroad On the specimen increases until a critical value, 
Ric Pomretcewea.smAct Chis point the crack rapidly propagates, 
anmomtt has been found that Kic Pomaemateriaal property, wsuen 
that in every structure when this condition is reached, the 
crack becomes unstable. 

MmvlLSec@eelascttc Materials, the crack is found to propa- 
feasemtor almoSt any load level. The velocities are infini- 
tesimally small for low loads and increaSe with increasing 
Toad= The question then becomes: at what velocity does 
mime rack propagate for any given load, not at what load 
level does rapid rade DiOwaceieton OCCUL. A stress inten- 
Pity factor POe@rene Viscoelastic case, cae can be obtained 
by the correspondence principle, where the stress intensity 
FACEOL; Ko FOGmiEMeweolastic case an the Laplace transform 
GCtiearn iseamnverted into the time domain to give the visco- 
Clacetemorreess Intensity factor, K, Sones Pm orina t Om 


relative to the specimen geometry, material, and boundary 


Goma tons < 
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eee ractire postulate Lor viscoelastic materials can 


be Stated in the £Lollowing way: 


- 


It 


Vv Vv 
a [K, (t)] (K, (t)] 


BOBBY 1 BODY? 2 


enn eva Rome) = [c (t)] 


BODY 1 BODY 2 


ime rommilReia 2) COoOnaucted constant stress and constant 
Mpeplacement tests on filled viscoelastic materials in an 
Mereem pe tO substantiate this relation. He developed frac- 
Mmimcmehanacterizations by plotting stress intensity factors 
Mersusscrack velocities for each test. His results showed 
mremconstant displacement and constant stress fracture 
emetacteriZzdt1ons diverging with increasing stress intensity 
mieror. ihis disparity was attributed to the dewetting 
memaviror of the material. 

Tim@emscecuayemakes a £Lurther attempt to prove the fracture 
postulate uSing an unfilled Visc@cleasctice material, since 
[Wome lrtetime hs eo NOE Experience dewetting. Fracture 
Si vacuerizgations were developed, using the same procedure 
as neediest s: from three experimental tests: ENG seOns Fant 
Pete acse weenie GCONStLant GAisplacement case, and the constant 
Serieement Late case. A cycling application of the con- 
stant displacement rate test was ee performed as an 
ag@edamuaistory Common to viscoelastic analysis and as a fur- 
Piemmcestmor the fracture postulate. These test results 


WetenwL nen Comeared, to evaluate the fracture postulate. 





’ rie Peo PRE MENTAL EOUILPMENT 


aie tool oPECIMENS 
The unfilled viscoelastic material used in testing was 
a carboxy terminated, polybutadiene (CTPB) polymer composed 
of 94.97% Butarez and 5.03% Hx868. The samples were fabri- 
foapeaubpy the Naval Weapons Center, China Lake, California, 
The specimens were prepared for testing by cutting a 
three-fourths inch horizontal crack in the sample as shown 
in Figure il. Emeemmens coulda then be used for crack lengths 
up to three inches with less than 1% difference between the 
far-field stress in the undisturbed region and the average 


mere aembacecadmom metmeross section [Ref. 1]. 


; | Wooo FRAME \ f , 
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a la Viscoetastic S7ATERIAL fZ|—— 
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Migtecml. — Geometry of Test Specimens 





Bi, LABORATORY HARDWARE 

Experimental testing was done with a table model 
Universal Instron. Pecemperature chamber was built to 
meeeene Instron so that specimen temperature could be 
megulated during the test (see Figure 2). The Instron 
Pee duipoed to directly conduct constant displacement, 


Mees cAnememoammacemMent rate, and cyclic tests. 





Figure 2. Universal Instron with Temperature Chamber. 


For the constant load test, the apparatus was adapted 
(Figure 3) so a constant load could be suspended beneath 
the temperature chamber. Since this case required a 
uniform stress state on the loaded boundary, the apparatus 


was built so that the center of gravity of the applied 


10 





load could be maintained under the center of the uncracked 
length of the specimen. This provided a constant strain 
along the length of the specimen, insuring a uniform stress 


boundary condition, 





Figure 3. Universal Instron and Temperature Chamber 
Adapted for Constant Load Test 

Specimen temperature was monitored with five thermis- 
tors each located one-eighth inch behind the sample and 
Spaced so as to provide a thorough coverage of the 
material surface. Their locations are shown in the photo- 
graph of Figure 4. 

Even though the temperature coils were exposed and 
located relatively close to the specimen ends, minimal 


Mempewdmure GradientS occurred in the area of the test 


Wal 










THERIMISTOR 
NUMG ER F. 


paresteteetete! 
W a 


Figure 4. Temperature Chamber. 


Sample. The uniform temperature in the horizontal direc- 
tion was attributed to the large metal fixture in the 
center of the chamber which acted as a heat Sink, provided 
much thermal inertia, and tended to stabilize the tempera- 
ture in its immediate vicinity. The temperature chamber 
had a substantial vertical temperature gradient, but this 
feature was not critical since the specimen dimension in 
this direction was small. 

Temperature chamber statistics are given in Table l. 


The "high" and "low" data represented the test specimen 


in 





environment at the maximum and minimum temperatures with- 
Pmcnien chamber during One cycle. PH 1.S high-low-high 
chamber cycle PiwIienhiverminutes for the 15, 30 and 
50 degree temperatures, while the compressor ran for a 
eager portion of the cycle as the temperature decreased. 
At 70 degrees, the chamber was basically maintaining room 
temperature but providing a chamber environment. At zero 
meme c>s the compressor was continually running providing 


pee latively constant temperature. 


THEA NISTOR NUNGER 


FPOOM TEMPERATURE 


fo | compaction cumrtnatey puma 
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tiaooe tee Lemperature Chamber Statistics 






















TEST TEMPERATURE, °F 
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- hi. Pee eee tA RESULTS 


A. SeONoTANT LOAD CASE 
The viscoelastic stress intensity factor (opening 
mode) for the constant load case was found by Lindsey 


meet. 1] to be: 


V _ Bit) IE 
co V Att aCe ee vr, 


Tpeomeeest COnsisted Of placing a step load on the 
sample, which was then maintained constant throughout 
meme Lest. Load levels required to produce crack propa- 
Goeren created specimen strain levels from 20% to 408%. 
To provide a uniform stress on the loaded specimen 
poaiemcary and tO avoid introducing a moment at the crack 
mods the rack propagated, the center of gravity of 
the applied load was continually moved so as always to 
be under the center of the uncracked length of the speci- 
men. Sbeweeciackevelocitics made this method tractable. 

SENecemtme specimen wasS finite in length, an increase 
miecktack Length resulted ina decrease in A(t) and an 
aoeseciated ancréease in material strain. This increased 
the sample width and, by Poisson's effect, reduced the 
i eomuameci ness. | With Increasing crack length, then, 
Giemotress intensity Lactor increased, resulting ina 


history Of Stress intensity factors for one applied load. 


14 





For the constant load case, the increase in crack 
mairoeclty with an inerease in crack length was small. 
Before, the crack velocity waS experimentally observed, 
pumamaveraged about four times, during each test. The 
stress intensity factor was computed based on the average 
mmek length for a particular crack velocity. Stress in- 
meemsity factors versus crack ve leek elas WeWe wtienm plotted 
merespecific testing temperatures. Since the visco- 
mitstic Material was thermoreologically simple [Ref. 2], 
each temperature curve was ultimately shifted into a 
miteotnh Master Eracture characterization curve. 

Individual temperature test results are given in 
meamre 5, and the Master fracture Characterization curve 
Menmiven in Figure G6. Test data had relatively little 
Seatter, and the experimental temperature shift factors 


micamgood COrrellation with relaxation modulus shift 


mactors [see Figure 6]. 


B. CONSTANT DISPLACEMENT CASE 
The viscoelastic stress intensity factor (opening 
fe@ey for a constant displacement case was found by 


Lindsey [Ref. JL || exes Jeter 


Vv 4 
K) (©) = Wapcerm Yo Prep!) ee 


To obtain a constant displacement boundary condition, 
the specimen boundary was displaced at a rate of two 


Micnes peLeminute Untill reaching the desired displacement, 
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ve" Then, by waiting ten times the loading time, the ramp 
tests Simulated step displacement loading for which the 
femess intensity factor reduced to Equation (2). As time 
mieteased after loading, the stress intensity factor and 
corresponding crack velocity decreased as a result of 
meet ion. This yielded a history of fracture data from 
eme test. 

meet the constant load case, the crack velocity was 
experimentally observed and averaged over each time incre- 
ment. Stress intensity factors were computed based on 
the average relaxation modulus for a particular crack 
Peewu@ecaty, Stress intensity factors versus crack velocities 
temenmplotted as a function of temperature, and finally 
Pmmweeed tO a mMaSter fracture characterization curve. 

fimaividual temperature test resultS aré given in 
Pane 7 and Figure 8, and the master fracture characteri- 
Moaemem curve iS given in Figure 9. Test data had a 
Deemerate degree of scatter, and experimental temperature 
Sreee@eetactors had some disagreement with relaxation 


modulus shift factors [see Figure 9]. 


cy GorscitaAanrt DISPLACEMENT RATE CASE 
teemvasecoclastic Stress intensity factor (opening 
mode) for a constant displacement rate case was found by 


Pamnasey (Ref. 1j] to be: 


ce E 
Ky (t) = (o> Rs, (1 S cio )e/|) ot E,t! (3) 


iis: 
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Figure 7. Fracture Characterization Curves for the 
Conctantepasplacement Case 
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Figure 8, Fracture Characterization Curves for the 
Constant Displacement Case 
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where EO, were determined by the procedure given in 


feneG SOlid Bropellant Mechanical Behavior Manual [Ref. 3]. 


# 


The constant displacement rate test utilized the 
moet ron in itsemost straight-forward application. One dis- 
placement rate was selected for the Instron crossarm, and 
this rate was maintained throughout the complete test. 

This test differed from the constant load test and the con- 
stant displacement test in two-.ways. First, one test 
generated a wide range of crack velocities, prohibiting 

mime averaging of crack Saar during Salve test < Second ly, 
this extensive history of data from one test resulted in 
@mly three different constant displacement rates at one 
maimeerature being required for sufficient fracture charac- 
jerilzation. 

Tomaccurately represent experimental data, crack 
lengths versus time were plotted. Practimesenaracteriza= 
tion was then done in two wayS. Pues teres aen velocities 
were graphically obtained by measuring the Slope of the 
@rack length ores time PlLOts ,,and coreesponding Stress 
Ln Wenn es factors were computed. This method permitted 
a eet Comparison with previous results: 

oem. recognizing the inaccuracy inherent in 
meaonmcalliy obtaining crack velocities, an incremental 
technique was employed. mhe crack Eee Suse ane joe eke 
was divided into discrete time increments. Beginning at 
Gnomieset time imerement, a stress intensity factor was 


Computed from Equation (3), and a corresponding crack 


Die, 





mebocity was Obtained from the master fracture characteri- 
moteswon curves. The current crack length was then obtained 
mrom the PEAK VelOCrtiy oy asSsumeng Jt to be ‘constant over 
the time increment. This process was continued and the 
resulting crack length-time history was compared with the 
experimental crack length-time history. 

Experimental results based on graphically obtaining 
meme rack velocity are given in Figure 10, labeled "Con- 
Stant Displacement Rate." Crack length versus time plots 
are given for the three constant displacement rate tests 
meamergure Jl, Figure 12, and Figure 13, and are labeled 
meee krLimental, Constant Displacement Rate.” A smooth, 


continuous crack length-versus-time curve resulted from 


the constant displacement rate test. 


Dr GreciiC CASE 
The basic equation for a viscoelastic stress intensity 
mieror (Opening mode) for the cyclic case was given by 


iamasey [Ref. 1] as: 


—— tc 
V 4 av 
K Ct oe —— f Oo 
- —— 4 
1 3bT - Ee sie) ae atm (4a) 
= ti) 

= = >, Ab 
a te BE. e i +E, (4b) 


temevcolic case was an extension of the constant dis- 
placement rate case as one crossarm rate was maintained 
mimedohettethemcemplete test. The direction of crossarm 


travel was reversed, however, when a selected displacement 
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Figure 13. Crack Length versus Time for the Constant 
Displacement Rate Case 
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was reached. This resulted in a linear, piecewise con- 


tinuous, even function which could be approximated by a 


fjmomrier series, 


oN 
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mere 14, Basic Cycle for the Cyclic Test. 
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Substituting Equation (ob) and? Equatwonms14)) nto 


Equation (4a) resulted in: 
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Equation (6) was solved on the computer using 14 terms in 
the Fourier series. 

Experimental data were plotted as crack length versus 
EIme . This was then compared with the constant load case 
and constant displacement case using the incremental pro- 
eedure described in the constant displacement rate case. 

perereck length-versus-time plot is given for a cycling 
meses in Figure 15 (also shown in Figure 16), labeled, 
Meeaeerimental, Cyclic Case." This test produced a stair- 


peep crack length history. 
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, TV DISCUSS LON SOF SRESULTS 


he constant load, constant displacement, constant dis- 
Merc cment rate, and cyclic tests were diverse in nature, 
each having several positive aspects and some having minor 
feeeectencies. Each of these four tests probed different 
aspects of material behavior, and collectively they provided 
@ecnorough and comprehensive examination of the fracture 
mectulate. 

Mie constant load case had two positive features. The 
meremoccurred over a relatively long period of time, simpli- 
Mac laboratory technigues and stabilizing test conditions, 
pmeetlne theoretical stress intensity factors were indepen- 
Gent ot external relaxation data. Detwmacting sfram tis 
test, however, was the mechanical apparatus uSed to apply 
aeeonsStant load. Unless carefully manipulated, it would 
iidw6ce a moment at the crack tip. 

che constant ac aaak case had several attractive 
Reacures . Stress intensity factors were calculated using 
mmr tlonwdata obtained, from constant displacement, tests. 
Sumi@e the test material was an elastomer, the reduction in 
Poem velocityadue to ,relaxationsoccurred over a substan- 
mEwulyebonger time than a similar test using solid propel- 
ete ees lemiewlete [Ref ,<.02)] LEoundethis test vigeual ly jampos- 
sible to use for propellants because rapid relaxation 


caused the crack velocities to go to zero before data 


a2 





could be taken. Finally, the Instron was capable of a 
Semsistent application of the constant displacement caSe. 
Prewconsctane displacement rate case utilized the 
MiseGon in Les most fundamental mode of operation. Since 

mewee.estS were required for fracture characterization, 
these tests were repeated. Excellent consistency was 
obtained between like tests. A primary disadvantage was 
Mime tnme Stress intensity factor calculation required 
meWaxation data. The only relaxation data available were 
mmosce Obtained from a constant Strain case. 

ive meOmmenes pertaining to the constant displacement 
memmeercase were relevant to the cyclic case. The cyclic 
case was a more thorough test, however, as disparities 
between experimental and theoretical results were accumu- 
Marea OVer G€ach cycle, and a discrepancy not apparent in 
one cycle (equivalent to one constant displacement rate 
Meat lewas apparent after several cycles. 

Hitec) er tne taster tracture characterization 
Sleves Given in Figure 10 graphically compares the con- 
Meaiteaisplacement, Constant load, and constant displace- 
ment rate cases. iiiseshows a correlation of results 
among the three sole mat sehr Om band Of +/% measured 
em the “vercical . Pigment Mowies l2, and Figure 13 com- 
pares the three constant displacement rate cases with the 
Pee boacdmwecace and the Constant displacement case 
Moangeehe Imeremental technique. These results were con- 


Sistent with Figure 10, as moderate correlation was 
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Soetained in all test cases. Page. 15> and Figure 16 com- 
Meme the cyclic case with the constant load case and the 
eonstant Pe Mee enone Casey cespectively, using the incre- 
Pameale technique. The cyclic test had good correlation 
Mmmeietme constant load test through four cycles. Le-would 
memarerficult tO maintain this agreement through a large 
number of cycles, however, as the incremental process 
mesultead in the error in one cycle being added to succeed- 
maa cycles. Stomianteeseat lon Of the cyclic test is 
mimeaect trated in Figure 16 where the constant displacement 
and constant displacement rate cases became further sepa- 
iteead with increasing cycles, because only a moderate 
agreement existed between the two test cases in any one 
evycle. 

Miiemanialyerecal theory used in obtaining viscoelastic 
stress intensity factors was based upon linear theory in 
that the change in the sample dimensions, when it was 
loaded, was assumed to be small. PimaecrlateeeesSt cases, 
however, Be rich s experienced from 20% to 40% strain, 
saueiae Semsomwdcaab le yariation im Sample width and thick- 
ness. This represented a significant departure from the 
linear case. 

A second departure from the analytical PiMmeOty Ceccus red 
wmecimeem@emeclack developed a finite crack tip radius. The 
Slasslealetieony for the stress intensity factor was 
flew leped £On an infinitely sharp crack. When the unfilled 


viscoelastic test specimens were loaded, however, the crack 
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Llp waS approximately elliptic. Figure 17 shows two 
Specimens Peano rerereiemserarms under typical loading 
Peewee Onweethese photographs indicate not only a finite 
emack tiperad@us butealso a@@echange in radius with a 
Sorresponding ehange in material strain. Since the con- 
aoe Ol aeStress I nNtenslttey stactLOr aS a Lunction of crack 
tip radius has not been developed, one can only specu- 


[Mace on its effect. 


oa) 








Figure 17. 


.275, p = .053, T = 70 


Typical Crack Geometries. 
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V. CONCLUSION 


wire Moderate Correlation among the constant load, con- 
Bewiewadlsplacement, constant displacement rate, and cyclic 
test cases tended to support the fracture postulate and 
gave credence to Hertzler's argument concerning dewetting. 
Roeeeseoulto, however, lacked absolute agreement as there 
Waoewa 15% Spread among the fracture characterizations for 
the four test cases. This disparity was primarily attri- 
buted to the nonlinearities inherent in the experimental 
meses. The effect of a finite crack tip radius on stress 
Mmimemstty Lactor waS not known, but its influence did not 
epeam EO be large since only a small range of crack tip 
radii were present in the experimental tests. 

The fracture postulate would appear to be even more 
applicable to propellant, assuming the dewetting phenomenon 
Geomra be Incorporated into the fracture theory, because 
the relatively Rite nature Cie sOPec Plant and Ene coOn— 
meeaimts imposed by the rocket motor structure would limit 
material ee within the linear assumptions used in 


ne -Eracture postulate. 


37 





REFERENCES 


Hawa Postgraduate Sehool = 57L1/73031A, Viscoelastic 
Pik si setietmedutor Rocket Motors Containing 
Gracksumoy, Gerald H- Lindsey, March 1973, 


Hertzler, Charles Miller, Viscoelastic Fracture 
Clasoerehimotteoneor a SOlid Propellant, A. E. 
Thesis, Naval Postgraduate School, Monterey, 


Coll homudied, sedldne 197 2 < 


Gnemmecal Propulsion Information Agency Publication 21, 
Bene Gwooome opellant Mechanical Behavior Manual, 
iiemiOonimiiepk imac Umiversity Applied Physics 
PioematOmy mol ilver spring, Maryland, September 
IMSS) S54 


a3 





, Thieves LR RUT PON LIST 
INIGr: 


Defense Documentation Center 
Gameron Station 
ele Sci tat ae lca nia 9 22314 


io rary pecode~O212 
Naval Postgraduate School 
Monterey, California 93940 


Associate Professor G. H. Lindsey 
Department of Aeronautics 

Navel Postgraduate School 
Monterey, California 93940 


Eeeutenant Robert Edwin Kapernick, USN 
meso Villa Circle 
Man Mina Calaforn ia 


Chairman 

Department of Aeronautics 
Naval Postgraduate School 
Monterey, California 93940 


39 


Copies 











SECURITY CLASSIFICATION OF THIS PAGE (Kien Data Entered) 


= te ow 


REPORT DOCUMENTATION PAGE 


REPORT NUMBER 2. GOVT ACCESSION NO. 
















READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


RECIPIENT'S CATALOG NUMBER 





4. TITLE (and Subtitie) 5. TYPE OF REPORT & PERIOD COVERED 
Master os. Thesis 


Eevarwation of Mew rrackuY Postulate 
: December 1973 


moemeyascoelastic Materials 
6. PERFORMING ORG, REPORT NUMBER 





7. AUTHOR(a) 8. CONTRACT OR GRANT HUMBER(e) 
Robert Edwin Kapernick 


PERFORMING ORGANIZATION NAME AND ADDRESS ] 10. PROGRAM ELEM 
AREA & WORK U 


ENT, PROJECT, TASY 
NIT NUMBERS 


Naval Postgraduate School 
memterey, California 93940 


CONTROLLING OFFICE NAME AND ADDRESS 12. REPORT DATE 
Deeemiber 1.973 


Naval Postgraduate School 
13. NUMBER OF PAGES 


momeerey, California 93940 









4] 

14. MONITORING AGENCY NAME & ADDRESS(I! diiferent from Conirolling Office) 15. SECURITY CLASS. (of thia report) 
Naval Postgraduate School Unclassified 
Monterey, California 93940 

(Sa, DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 





16. DISTRIBUTION STATEMENT (of thie Report) 


Paeezoved £Or public release; distribution unlimited. 


17. DISTRIBUTION STATEMENT (of the ebhetract entered in Black 20, if different from Report) 





16. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue on reveree olde if neceeeary end identify by block number) 


Eracture 
wrseoelastic Fracture 
meilastomer Fracture 





20. ABSTRACT (Continue on reveroo aide if neceeeery end identify by block number) 


Bxperimental studies were conducted on unfilled viscoelastic 
Materials to evaluate the fracture postulate which states that 
Pqermevelociny 15 4a unique function of stress intensity 
hae LOL. Romcomtiits,s Fracture Characterization tests were con- 
ieee Cureoumconotant Load, constant displacement, and constant 
lestecnenmtemdte histories. A cyclic study, a more sophisti- 


Sect pleat lonmor the constant displacement rate case, was 
0 eS EE Pe SB ES BA TT I 





EL EE 0 REO EO A ES Ee FI 
DD pee; 1473  evition oF 1 Nov 65 1S OWSOLETE 


o S/N 0102-014- 6601 —<—$—$—$— rt De Ba 
(Page 1) : : SECUTITY CLASSIFICATION OF THIS PAGE (Whon Deta Entored} 


40 








SECURITY CLASSIFICATION OF THIS PAGE(When Deta Entered) 








meso conducted as a further evaluation of the fracture postulate 

Micmiialy emo MNecory Wsed 1n Obtaining viscoelastic stress 
intensity factors was based upon linear theory and sharp crack 
geometries. Moderate correlation was obtained among all test 
Sises, -NOLWLENStanding the presence of Significant specimen 
Peiatinoeand “Linlte Crack tip radii. 





DD ; Form. 1473 (BACK) 
an FE ee or ene 
S/N 0102-014-6601 | SECURITY CLASSIFICATION OF THIS PAGE(When Deta Entered) 
4 1. 




















